Accumulation of three-body resonances above two-body thresholds 
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We calculate resonances in three-body systems with attractive Coulomb potentials by solving the homogeneous 
Faddeev-Merkuriev integral equations for complex energies. The equations are solved by using the Coulomb- 
Sturmian separable expansion approach. This approach provides an exact treatment of the threshold behavior 
of the three-body Coulombic systems. We considered the negative positronium ion and, besides locating all 
the previously know S-wave resonances, we found a whole bunch of new resonances accumulated just slightly 
above the two-body thresholds. The way they accumulate indicates that probably there are infinitely many 
resonances just above the two-body thresholds, and this might be a general property of three-body systems with 
attractive Coulomb potentials. 

PACS numbers: 34.10,+x, 31.15.-p, 02.30.Rz 



The most common method for calculating resonant states 
in quantum mechanical systems is the one based on the com- 
plex rotation of coordinates. The complex rotation turns the 
resonant behavior of the wave function into a bound-state- 
like asymptotic behavior. Then, standard bound-state meth- 
ods become applicable also for calculating resonances. The 
complex rotation of the coordinates does not change the dis- 
crete spectrum, the branch cut, which corresponds to scatter- 
ing states, however, is rotated down onto the complex energy 
plane, and as a consequence, resonant states from the unphys- 
ical sheet become accessible. By changing the rotation an- 
gle the points corresponding to the continuum move, while 
those corresponding to discrete states, like bound and reso- 
nant states, stay. This way one can determine resonance pa- 
rameters. In three-body systems there are several branch cuts 
associated with two-body thresholds. 

In practice, the complex rotational technique is combined 
with some variational approach. This results in a discretiza- 
tion of the rotated continuum. The points of the discretized 
continuum scatter around the rotated-down straight line. So, 
especially around thresholds it is not easy to decide whether 
a point is a resonance point or it belongs to the rotated con- 
tinuum. Moreover, variational methods approach states from 
above, so resonances slightly above the thresholds may easily 
get lost. 

Recently, we have developed a method for calculating res- 
onances in three-body Coulombic systems by solving homo- 
geneous Faddeev-Merkuriev integral equations 1 1 ] using the 
Coulomb-Sturmian separable expansion approach |2]. As a 
test case, we calculated the resonances of the negative positro- 
nium ion. This system has been extensively studied in the past 
two decades and thus serves as test example for new methods. 
We found all the 12 S-wave resonances presented in Ref. |3] 
and observed good agreements in all cases. 

We also observed that in case of attractive Coulomb interac- 
tions the Faddeev-Merkuriev integral equations may produce 
spurious resonances |4], which are related to the somewhat 
arbitrary splitting of the potential in the three-body configu- 



ration space into short-range and long-range terms. We could 
single them out by changing those parameters. We succeeded 
in locating 10 more resonances in the same energy region, all 
of them are very close to the thresholds. These new reso- 
nances were published in Ref. U 

As our skill in applying our method developed we located 
more an more new resonances just slightly above the two- 
body thresholds. They are all aligned along a line in the com- 
plex energy plane pointing toward the thresholds. It seems 
that there are infinitely many resonances accumulating at the 
two-body thresholds. Since our method is relatively new we 
briefly outline the basic concepts and the numerical tech- 
niques, specialized to the e - e~e + system (further details are 
in Refs. BH) 

The Hamiltonian of a three-body atomic system is given by 

H = H° +v? + v% + v%, (1) 

where H° is the three-body kinetic energy operator and 
denotes the Coulomb potential in the subsystem a, with a = 
1,2,3. We use throughout the usual configuration-space Ja- 
cobi coordinates x a and y a , where x a is the coordinate of the 
(/3, 7) pair and y a connects the center of mass of (/?, 7) to 
the particle a, respectively. Thus v^, the potential between 
particles j3 and 7, depends on x a . 

The Hamiltonian Q is defined in the three-body Hilbert 
space. The three-body kinetic energy, when the center-of- 
mass motion is separated, is given by 

H° = h° Xa + h° ya = h° X0 + h% 3 = h x ^ + h ^ , (2) 

where h° is the two-body kinetic energy. The two-body poten- 
tial operators are formally embedded in the three-body Hilbert 
space v c = v c (x)l y , where l y is a unit operator in the two- 
body Hilbert space associated with the y coordinate. 

In Merkuriev's approach to the three-body Coulomb prob- 
lem 1 1 ] the Coulomb interaction is split, in the three-body con- 
figuration space, into short- and long-range terms 
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FIG. 1: i/ s ' and i/ s ' for an attractive Coulomb potential. 

where the short- and long-range parts are defined via a split- 
ting function: 

«i s) = v«((x a ,y a ) (4) 

«i = «a [W(*a,lte)]- (5) 

The splitting function £ is defined such that 

i im C(*,y) = {l' if N<-o(i + MM)^, (6) 

x, y ^oo sv ,y; \ 0, otherwise, 

where £0,2/0 > and ^ > 2. So, in the region of three- 
body configuration space where particles (3 and 7 are close to 



each other Va 



and Va ~ 0, otherwise viz ~ w£ and 



0. Usually the functional form 

C(z, y) = 2/ {1 + exp [{x/x ) v /(l + y/y )}} , (7) 

is used. Typical picture for and i>(" are seen in Fig.^ 

In atomic three-particle systems the sign of the charge is 
always identical for two particles. Let us denote those two 
particles by 1 and 2, and the third one by 3. In this case v% 
is a repulsive Coulomb potential which does not support two- 
body bound states. Therefore the entire can be considered 
as long-range potential and the Hamiltonian can formally be 
written in a form which looks like an usual three-body Hamil- 
tonian with two short-range potentials 



H = ffW 



,( s ) _lJ s ) 



where the long-range Hamiltonian is defined as 



(8) 



(9) 



Then, the Faddeev method is applicable and, in this particu- 
lar case, results in a splitting of the wave function into two 
components 



(10) 



The components are defined by \ip a ) = G^(z)v a 
where a = 1, 2 and G^> (z) = (z — H^')~ x , z is a complex 
number. 

In the cases of bound and resonant states the wave-function 
components satisfy the homogeneous two-component 
Faddeev-Merkuriev integral equations 



|V>t) = G«(*>W|V> 2 > 
|^ 2 ) = G®(z)vP\ih) 



(11) 
(12) 



at real and complex energies, respectively. Here Ga is the 
resolvent of the channel long-ranged Hamiltonian Go? (z) — 
- h2 ] )-\ where H ( J ] = ffW + 

Further simplification can be achieved if we take into ac- 
count that particles 1 and 2 are identical and indistinguishable. 
Then, the Faddeev components and ^2), in their own 
natural Jacobi coordinates, have the same functional forms 
{xiyi\ipi} = (x2j/2|V , 2)- On the other hand | ^2) = pV^i), 
where V is the operator for the permutation of indexes 1 and 
2 and p = ±1 denotes the eigenvalue of V. Therefore we can 
determine from the first equation only 



|V-'i) = Gfv^pP^i). 



(13) 



It should be noted, that so far we did not make any approxi- 
mation, and although this integral equation has only one com- 
ponent, yet it gives full account both of asymptotic and sym- 
metry properties of the system. 

We solve Eq. Jl 3i by using the Coulomb-Sturmian separa- 
ble expansion approach. The Coulomb-Sturmian (CS) func- 
tions are defined by 



il) = 



(n + 2l + iy. 



1/2 



(26r) i+i cxp(-br)Lf!+ L (2br), 

(14) 

n and I are the radial and orbital angular momentum quantum 
numbers, respectively, and b is the size parameter of the basis. 
The CS functions form a biorthonormal discrete basis 

in the radial two-body Hilbert space; the biorthogonal partner 
is defined by (r\nl) = (r\nl)/r. Since the three-body Hilbert 
space is a direct product of two-body Hilbert spaces an appro- 
priate basis can be defined as the angular momentum coupled 
direct product of the two-body bases 

\nvlX)i = \nl)i (g) \u\)i, (n, v = 0, 1, 2, . . .), (15) 

where \nl)\ and \vX)\ are associated with the coordinates x\ 
and yi, respectively. With this basis the completeness relation 
takes the form 



JY 



lim >^ \nvlX)-\ \(nvlX\ = lim if. 



71, U — 



N — >oc 



(16) 



Similar bases can be constructed for fragmentations 2 and 3 
as well. 

We make the following approximation on Eq. (1131 



\iP 1 )=G[ l h?v[ s) pVl?\i> 1 ), 



N\ 



(s) 

i.e. the operator v\ pV in the three-body Hilbert space is ap- 
proximated by a separable form, viz. 



(17) 



v\ s 'pV = lim If vf'pVl? « If v\ s >pVl? 

N^oo 
N 



\nvlX)i v[ s) i(n'u'l'X'\, (18) 



3 



where uj^ = i(ni>l\\v[ sS> pP\n' V V A')i. Utilizing the proper- 
ties of the exchange operator V these matrix elements can be 
written in the form — p x (— ) l i (nvl\\vi \n'v'l'\')2- 

With this approximation, solving Eq. (II 3I > turns into solving 
the matrix equation 



(19) 



-.(0 



for the component vector tp = i(nvl\\ip±), where G\ ' 

A unique solution exists if and only 



if 



D(z) = det{[Gi l) (z)]- L -v\ s> } = 



(20) 



So, to calculate resonances, we need to search for the zeros of 
determinant D{z) on the complex energy plane. 

The Green's operator G_i is related to the Hamiltonian 
Hi, which is still a three-body Coulomb Hamiltonian and 
seems to be as complicated as H itself. However this is not 
the case. The only possible two-body asymptotic configura- 
tion for Hi is when particles 2 and 3 form a bound states and 
particle 1 is at infinity. The corresponding asymptotic Hamil- 
tonian is 



Hi =H° 



(21) 



Therefore, in the spirit of the three-potential formalism 0], 
G_i can be linked to the matrix elements of G\{z) = (z — 
Hi) -1 via solution of a Lippmann-Schwinger equation, 



u 



■1) 



(22) 



(0 



where = i{nvl\\Gi\n' v'V \')i andf/j = i(nvl\\(v. 
v$)\n'i/V\') x . 

Now, what is remained is the calculation of the matrix el- 
ements G_i, since the potential matrix elements v_i and U_i 
can always be evaluated numerically. The Green's operator 
Gi is a resolvent of the sum of two commuting Hamiltonians, 



Hi = h x 



'in 



where h Xl = h Xi 



-wf {xi) mdh Vl = h yi , 



which act in different two-body Hilbert spaces. Thus, G_i can 
be given by a convolution integral of two-body Green's matri- 
ces, i.e. 



1 



= — ^ dz'g^(z-z')g^(z'), 



(23) 



where g Xl (z) = (z - h Xl )~ x and g Vl (z) = (z - h yi )~ x . 
The contour G should be taken counterclockwise around the 
continuous spectrum of h yi such a way that g Xl is analytic 
on the domain encircled by G. With the contour on Fig.[2]this 
mathematical condition is met even for resonant-state energies 
with z = E — iT/2. The corresponding CS matrix elements 
of the two-body Green's operators in the integrand are known 
exactly and analytically for all complex energies (see [[5[ and 
references therein). From this follows that all the thresholds, 




g x (E-iz-z') 



FIG. 2: Analytic structure of g xi (E + ie — z') g yi (z') as a function 
of z' , e = — r/2. The Green's operator g yi (z') has a branch-cut 
on the [0, oo) interval, while g xi (E + ie — z') has a branch-cut on 
the (— oo, E + ie] interval and infinitely many poles accumulated at 
E + ie (denoted by dots). The contour C encircles the branch-cut 
of g yi such that a part of it goes on the unphysical Riemann-sheet 
of g yi (drawn by broken line) and the other part detoured away from 
the cut. The branch-cut and some poles of g xi (denoted by full dots) 
are lying on the physical Riemann-sheet, some other poles (denoted 
by empty dots) are lying on the un-physical Riemann-sheet of g yi , 
respectively. Yet, the contour avoids the singularities of g xi . 



which correspond to the poles of g , are at the right location, 
irrespective of the rank N used in the separable expansion. 

To calculate resonances we have to find the complex ze- 
ros of the Fredholm determinant D(z) of Eq. J20l >. Between 



thresholds D(z) is analytic, therefore, due to a theorem of ho- 
momorphic functions 



uyt 
16] 



1 

2ttI 



D'(z)/D(z)dz = N c > 



(24) 



c 



where Nc is the number of zeros inside the contour C . By 
calculating J24i numerically we can decide whether a domain 
contains a resonance or not. 

We considered the S-wave resonances of the e~e~e + sys- 
tem. The resonances, found at the vicinity of thresholds, are 
seen in Fig. [5] The calculations were performed with three 
entirely different sets of parameters: xq — 18 and yo = 50, 
xq = 25 and yo = 50, xq = 5 and yo = 1000, while v = 2.1 
in all cases (the lengths are given in gq units). We found, that 
the results at N — 20 CS basis states and angular momen- 
tum channels up to I = A = 10 are well converged and they 
are rather insensitive for the choice of CS parameter b over 
a broad interval. The resonances displayed in Fig. [3] are sta- 
ble against the change of xq and yo parameters, they exhibit a 
remarkable 5 — 6 digits stability. 

We can see that the resonances are aligned along a line 
pointing exactly to the two-body thresholds. As we stretched 
the code and went closer and closer to the threshold we dis- 
covered more and more resonances. All of them were along 
the line. This indicates that the two-body threshold is an ac- 
cumulation point of the resonances, and probably there are 
infinitely many there. 

This conclusion is supported by our previous study of the 
e + + H system 01, where violent oscillations of the cross sec- 
tions just above two-body thresholds were found. Preliminary 
resonance calculations with the present method show that in 
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FIG. 3: Accumulation of resonances above the two-body thresholds. 



the e + + H system, where the violent oscillations were found, 
there are also accumulation of resonances. 
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